The number of open paths in an oriented 
^-percolation model 

Francis Comets*'^ Serguei Popov^'^ 
Marina Vachkovskaia*'^ 

February 1, 2008 



""^Universite Paris 7, UFR de Mathematiques, case 7012, 2, place Jussieu, F-75251 Paris 
Cedex 05, France 

e-mail: comets@math.jussieu.fr, url: http : //www. proba. jussieu. fr/^comets 

^Instituto de Matematica e Estatistica, Universidade de Sao Paulo, rua do Matao 1010, 
CEP 05508-090, Sao Paulo SP, Brasil 

e-mail: popovSime.usp.br, uri: http://www.ime.usp.br/~popov 

^Departamento de Estatistica, Instituto de Matematica, Estatistica e Computagao Cienti- 
fica, Universidade de Campinas, Caixa Postal 6065, CEP 13083-970, Campinas SP, Brasil 
e-mail: marinav@ime.unicamp.br, url: http://www.ime.unicamp.br/~marinav 

Abstract 

We study the asymptotic properties of the number of open paths 
of length n in an oriented £»-percolation model. We show that this 
number is e""'-^^^^^"'-^-*^ as n — > oo. The exponent a is deterministic, 
it can be expressed in terms of the free energy of a polymer model, 
and it can be explicitely computed in some range of the parameters. 
Moreover, in a restricted range of the parameters, we even show that 
the number of such paths is n -i/2^e"°(f)(l -ho(l)) for some nonde- 
generate random variable W. We build on connections with the model 
of directed polymers in random environment, and we use techniques 
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and results developed in this context. 
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1 Introduction and results 
1.1 Introduction 

In this paper we study the number of open paths in an oriented ^^-percolation 
model in dimension 1 + d, or, equivalenly, the number of g-open path in 
an oriented percolation model. Consider the graph Z+ x Z*^, with Z+ = 
{0,1,2,3,...}, and fix some parameter p G (0,1). To each site of this 
graph except the origin, assign a variable taking value 1 with probabil- 
ity p and with probability 1 — p, independently of the other sites. An 
oriented (sometimes also called semi- oriented) path of length n is a sequence 
(0, xo), (1, xi), (2, X2), . . . , (n, Xn), where xq = and are neighbours 

in Z"^, z = 0, . . . ,n — 1: viewing the first coordinate as time, one can think 
of such path as a path of the d- dimensional simple random walk. Fix an- 
other parameter g G [0,1]; the concept of ^^-percolation was introduced by 
Menshikov and Zuev in [T^ , as the occurence of an infinite length path with 
asymptotic density of Is larger of equal to g. As in classical percolation, this 
probability of this event is subject to a dychotomy [19] according to p larger 
or smaller than some critical threshold, which was later studied by Kesten 
and Su [15] in the asymptotics of large dimension. 

In the present paper, we discuss paths of finite length n, in the limit 
n ^ 00. An oriented path of length n is called g-open, if the proportion 
of Is in it is at least g. From standard percolation theory it is known that 
for large p there are 1-open oriented paths with nonvanishing probability, 
and from [19] that for any p one can find g larger than p such that, almost 
surely, there are ^open oriented paths for large n. However, the question of 
how many such paths of length n can be found in a typical situation, was 
still unaddressed in the literature. When finishing this manuscript, we have 
learned of the related work [T6] . 

In this paper, we prove that the number of different g-open paths of 
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length n behaves hke e^°'^^^^^~^°^^^\ where the exponent a{g) is deterministic 
and, of course, also depends on p and d. We prove that the function «(•) 
is the negative convex conjugate of the free energy of directed polymers in 
random environment. This model has attracted a lot of interest in recent 
years, leading to a better - although very incomplete - understanding. We 
will extensively use the current knowledge of thermodynamics of the poljTiier 
model, and the reader is refered to [5| for a recent survey. This will allow 
us to obtain, when d > 3, the explicit expression for a{g) in a certain range 
of values for g depending on the parameters p and d. The reason for this 
remarkable fact is the existence of the so-called "weak-disorder region" in the 
polymer model, discovered in [13] and [2] : this reflects here into a parameter 
region where the number of paths is of the same order as its expected value. 

At this point the reader may be tempted to use first and second moment 
methods to estimate the number of paths. The first moment is easily com- 
puted, and serves as an upper bound in complete generality. The second 
moment is more difficult to analyse. However, it can be checked that in large 
dimension and for density close to the parameter p of the Bernoulli, the ra- 
tio second-to-first-squared remains bounded in the limit of an infinitely long 
path. This means that, under these circumstances, the upper bound gives 
the right order of magnitude with a positive probability. However, (i) this 
method does not tell us anything on a for general parameters, (ii) it fails 
to keep track of the correlation between counts for different values of the 
density. 

Our strategy will be quite different. We will study the moment generating 
function of the number of paths, which is not surprising in such a combina- 
torial problem. The point is that the moment generating function is simply 
the partition function of the directed polymer in random environment. This 
is a well-known object in statistical physics, its logarithmic asymptotics is 
well studied, and is given by the free energy. From the existence and known 
properties of the free energy, we will derive the existence of a and its expres- 
sion in thermodynamics terms. In the course of our analysis we will prove 
that the free energy, a convex function of the inverse temperature, is in fact 
strictly convex. This property is new and interesting for the polymer model. 

Moreover, in a more restricted range of values for g, we even obtain an 
equivalent for the number of paths which achieves exactly a given density 
of Is. This is clearly a very sharp estimate, that we obtain by using the 
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power of complex analysis, and convergence of the renormalized moment 
generating function in the sense of analytic functions. Certainly a naive 
moments method cannot lead to such an equivalent. 

1.2 Notations and results 

Now, let us define the model formally. Let r]{t,x),t = 1,2, ... ,x e Z*^ be a 
sequence of independent identically distributed Bernoulli random variables, 
with common parameter p G (0, 1), F{ri{t,x) = 1) = p = 1 — F{r]{t,x) = 0). 
We denote by {Q, A, P) the probability space where this sequence is defined. 
The vertex {t,x) is open if r]{t,x) = 1 and closed in the opposite case. A 
nearest neighbour path S in Z*^ of length n {1 < n < oo) is a sequence 
S ={St;t = 0,.. .,n), St e So = 0, \\St - St-i\\i = 1 for t = 1, . . . , n. We 
denote by Vn the set of such paths S, and by Voo the set of infinite length 
nearest neighbour paths. For S EVn, let 

n 

H^{S) = J2v{t,St) (1) 

i=l 

be the number of open vertices along the path S. 

In oriented percolation, one is concerned with the event that there exists 
an infinite open path S, i.e. 

Perc = {there exists S G Voo '■ vi^^ ^t) = 1 for all t > l} . 

It is well known [SI [12] that there exists Pc{d) G (0,1), called the critical 
percolation threshold, such that 

P(Perc)|>° ^! (2) 
^ ^ 1^ = if p < Pc{d). ^ ' 

For Q G {p, 1], Menshikov and Zuev [19] introduced ^^-percolation as the event 
that there exists an infinite path S with asymptotic proportion at least q of 
open sites, 

;0-Perc = (there exists S G Voo '■ liminf if„(5')/n > q\ . 

They showed that there also exists a threshold Pc{q, d) such that ([2]) holds 
with Q-'Perc instead of Perc (with the probability of g-'Perc being equal to 
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1 when p > Pc{Q,d)). Very little has been proved for ^-percolation. The 
asymptotics of Pc{g, d) for large d are obtained in [15] at first order, showing 
that d^^^pc{g, d) has a limit as d — >■ cxd, and that the limit is different from 
the analogous quantity for d-aiy trees. As mentioned in this reference, the 
equality pc{l,d) = Pc{d) follows from Theorem 5 of |18]. 

In this paper we are interested in the number of oriented paths of length n 
which have exactly k open vertices [k G {0, . . . , n}), 

QM = Card {SeVn:H^{S) = k} 

(Card A denotes the cardinality of A) and the related quantity given, for 
ge [0,1], by 

( Ca.Td{S eVn:H^{S)>ng} , g > p, 

Rn{g) = { , 

y Card {S eVn- H^iS) < ng} , g < p. 

Note that Qn{k), Rn{g) are random variables, that Rn{g) = J2k>ngQn{k) 
when g>p, and that Perc = n„{QnH > 1} = OniRni^) > 

In this paper we relate these quantities to the model of directed polymers 
in random environment. Central in this model is the (unnormalized) partition 
function Z„ = Zn{(3,f]) at inverse temperature /3 G M in the environment rj 
given by 

By subadditive arguments one can prove that 

(p{(3) = lim -ElnZ„ (3) 

n^oo n 

exists in M (E is the expectation under P), and by concentration arguments, 
that the event flo{P) defined by 

fio(/?) = I lim -lnZ„ = (4) 

has full measure, P(f2o(/?)) = 1, see e.g. |4j. The function ip is called the 
free energy, it is a non-decreasing and convex function of f3. Its Legendre 
conjugate 

ip*ig) = snp{Pg-ip{P);PeR}, (5) 
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is a convex, lower semi-continuous function from [0, 1] to M U {+00}, such 
that f*{g) > f*{p) = — ln(2(i) (indeed, f'{0) = p, as it will be shown 
later). Legendre convex duality is better understood by taking a glance at 
the graphical construction, e.g. figures 2.2.1 and 2.2.2 in here, on Figured] 
we illustrate how the functions and cp* typically look in our situation. The 
existence of the so-called time constants, 

+ r Hn{S) _ . H„{S) 
Q = hm max , g = hm mm , Ir-a.s. 

n^oo SeV„ n n^oo SeV„ Tl 

can be obtained by specifying a direction for the ending point Sn, which 
allows using subadditive arguments |l4j , and then summing over the possible 
directions. However we give here a short proof in the spirit of this work. Since 

exp{j3maxseVn Hn{S)} < Zn < (2(i)" exp{/?max5ep„ Hn{S)} , we have 

— - lnZ„ — - ln(2d) < max — — — ^ < — - In Z„ . 

nf3 13 ^ ' - sev„ n - np 

Taking the limits n —>■ 00 and then (3 — > +00, we see that is well-defined, 
and is in fact equal to the slope lim^^+ooV^(/5)//5 of the asymptotic direc- 
tion of if at +00. From standard properties of convex duality, the range 
of the derivative {d/ d[3){l/n) lnZ„(/5) a.s. converges to [g^ , g'^] in Hausdorff 
distance, and v'*(f?) < +00 if and only ii g E [g^ , g^]. For such g, we have 
ip\g) < 0. 

Theorem 1.1 For all g G [0, 1] with g ^ g^, g~, the following limit 

a{g) = lim — lni?„(^) (6) 
exists F-a.s. (possibly assuming the value —{00)), and is given by 

aig) = -^*ig) ■ 

Clearly, a is concave, with values in [0, ln(2(i)] U {—00} and a{p) = ln{2d). 
Note that, for all g & {g~ , g~^) and almost every 77, 

Rn{g) = expn[a(^) + o(l)] , as ri ^ 00 . 
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Figure 1: The function </? and its Legendre transform (/?* 
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Remark 1.2 

By convexity the function a is continuous on For now, it is not clear 

to us whether the hmit a^Q^—) should be equal to in the case p < Pc{d). 
In the case p > Pc{d), it is possible to show by subadditive arguments that, 
conditionally on percolation, the limit a(l) in ([6]) exists and is positive, but 
it is not clear to us whether a is continuous at 1. 

Let 

A(/3) = lnEe^''(*'") = In [l + p(e'^ - 1)] , \{p) = \{p) + \n{2d) , 

then KZn = exp{nX{(3)}. A direct computation shows that the Legendre 
conjugate A*(^) = sup{Pg — A(/3); /? G M} of A is equal to 

A*(^) = -\n{2d) + 0\n^ + {l-0)\nl^ (7) 

p 1 — p 

= gin — ^ — h (1 — f?) In ■ ^ 



2dp ' ' 2d{l-p)' 

The function (—A*) is important for understanding the rate a. (We recall that 
both functions depend on p, but we don't write explicitely the dependence.) 
Note that these two functions coincide at ^ = p and take the value ln{2d). 

Theorem 1.3 Let p e (0,1). 

1. We have the annealed bound.' For all g, 

a(f?) < -A*(f?) . (8) 

2. The function Q{g) + \*{g) is nonincreasing for g E [p, g^) and is non- 
decreasing for g G {g^,p]- 

3. The set 

Vip) = {geiO,l):aig) = -X*ig)} (9) 

is an interval containing p (here, "interval" is understood in broad 
sense, i.e., it can reduce to the single point {p}). 

4. In dimension d = 1, V{p) = {p}, i.e. the inequality in ^ is strict for 
all g ^ p. 
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V{p) 



Figure 2: Typical behaviour of the function a when d > 3 

5. In dimension d>3, V{p) contains a neighborhood of p. 

6. Let d > 3, and iTd be the probability for the d- dimensional simple ran- 
dom walk to ever return to the starting point. When p > TCa, then 
[p, 1) C V{p), so that the equality holds in ^ for all g G [p, 1). Sim- 
ilarly, when p < 1 — VTrf, then {0,p] C V{p), so that the equality holds 
for all Q G (0,p]. 

7. In dimension d > 2, if p < {l/2d), then supV(p) < 1. Similarly, if 
p> 1 - {l/2d), we have inf V(p) > 0. 

Remark 1.4 

(i) The annealed bound comes from the first-moment method, and most of 
the results stating that the equality a{g) = —\*{g) holds, are derived from 
the second-moment method. 

(ii) By transience of the random walk in dimension c? > 3, we have vr^ < 1. 
In fact, vTs = 0.3404... > 7r4 > 115... [23l page 103]. In particular, for 
p G {iTd, 1 — TCd), we have a{g) = —\*{g) for all g G (0, 1). 

The following property of the free energy of the directed polymer is 
interesting and seems to be new. 
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Theorem 1.5 The function (p is strictly convex on M, and the functions (f* 
and a are difjerentiable in the interior of their domains. 

We will obtain much sharper results for large dimension and g^s not too 
far from p. The reason is that the partition function Z„ behaves smoothly 
as n y oo. The almost-sure limit 

WooiP) = lim Z„(/5)e-"^(^) 

exists for all P, since the sequence is a positive (^„)„-martingale, where 
Gn = < n,x E Z*^}. So, let us now concentrate on the case 

of large dimension, d > 3. When j3 belongs to some neighborhood of the 
origin (known as the weak disorder region), the limit Woo is strictly positive 
a.s. In a smaller neighborhood of the origin, the limit can be expressed as 
a (random) perturbation series in 

£2 [22]. Moreover, the convergence holds 
in much stronger sense, namely, in the sense of analytic functions [7]. We 
will use strong tools from complex analysis, as it is classically done to obtain 
limit theorems for sums of random variables [20] . 

Theorem 1.6 Assume d > 3. There exist a neighborhood of p m M 
and an event ^2 with full probability such that for every sequence kn with 
kn/n g E Us and all rj e fl2, 

where o(l) tends to as n oo, and f3{g) = In ^^j-^- The neighborhood Us 
is contained in V{p), hence we have a = —A* with A* given by ([^. 

We note that the leading order is deterministic, but the prefactor is random 
(as Woo), depending on the particular realization of the Bernoulli field. This 
theorem is a corollary of a more refined result (Theorem 13.21) . which can be 
found in Section [3l This will be proved by complex analysis arguments, con- 
sidering the Fourier transform of Hn under some (polymer) measure. Fourier 
methods are quite strong, they are used in a different spirit in [Ij to obtain 
sharp results on the polymer path itself for small (3. The disadvantage is that 
we have to restrict the parameter domain. It would be tempting to use only 
real variable techniques as in the Ornstein-Zernike theory for the Bernoulli 
bond percolation [3], but we take another, shorter route. 



10 



Remark 1.7 

The model is also interesting with real-valued ri{t,x) with general distri- 
bution. This is motivated by first-passage time percolation. Our results 
at the exponential order remain valid for variables with exponential mo- 
ments. For the case of the Gaussian law, we mention the recent preprint 
[T7] on the so-called REM conjecture: it is proved that the local statistics of 
{Hn{S)] S G Vn) approach that of a Poisson point process, provided that one 
focuses on values distant from the mean Ei?„ by at most o{n^''^). 

We can interpret our last result in this spirit. In our case, (if„(5'); 5* G Vn) 
spreads on the lattice, and natural local statistics of the energy levels are the 
ratios Qnikn) /^Qn{kn) ■ For d>3 and kn ^ ng E U3, 

Qn{kn)/EQn{kn) ~ W^i^ig)) 

since KWooiP) = 1- W^e emphasize that here the energy level /c„ is of order n 
(far from the bulk), and that the limit is not universal but depends on the 
lattice and the law of the environment rj. 



2 Logarithmic asymptotics 

2.1 Proof of Theorem [TTI 

We start by introducing some probability measures. 

Let P be the law of the simple random walk on Z'^ starting from 0, i.e. the 
probability measure on the space Poo of infinite paths making the increments 
St — St-^i independent and uniformly distributed on the set of 2d neighbors 
of G Z'^. Observe that the restriction of P to paths of length n is the 
normalized counting measure on Vn, and so the partition function takes now 
the familiar form {Ep is the expectation with respect to P) 

Zn = {2drEp[exp{f5Hn{S)}]. (10) 

The law z/„ = z/^ of (l/n)if„ under P, given by ^'n({f?}) := PiHn{S) = ng), 
is such that 

^n(M) = ^||P^ if nge{0,l,...,n}. (11) 
We extend to a probability measure on M that we still denote by z/^, 

^niA) = E,eA,n,e{0,l,...,n}'^n({rf), A C M. 
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All what we need to obtain the proof of Theorem 11.11 is to prove that 
Vn obeys an almost sure large deviation principle, see Proposition 12. II below. 
Recall first the event fio(/3) from (jlj), and define the event = H/jgq ^o(/3). 
Then, we have P(fio) = 1? and on this event the convergence (jlj) holds for 
any real number (3 by convexity and monotonicity. 

Proposition 2.1 The function 

I{q) = \xi{2d) + i^*{Q) e [0, ln(2d)] U {+00} 

is lower semi- continuous and convex on [0, 1]. Moreover, for all rj & Qq the 
sequence (z/„,n > 1) obeys a large deviation principle with rate function I. 
That is, 

(i) for any closed F C [0, 1], we have 

lim sup In z/„ (F) < — inf I{g), 

n— >oo Q&F 

(a) for any open (in the induced topology on [0, 1]^ G C [0, 1], we have 

liminf In z/„(G) > — inf I{g). 

Now, we first finish the proof of Theorem 11.11 and then prove the above 
proposition. 

Proof of Theorem \l.l\ Assume that q G [p, q'^) and r] G VLq. Applying (i) of 
Proposition 12.11 with F = [g, 1] and using ffTTj) together with the fact that 
Q > p, we see that the limit in ([6]) is not larger than ln{2d) — I{g) = a{g). 
Applying (ii) of Proposition 12.11 with G = {g + e,l] {e > 0) and using the 
fact that g > p, we see that the limit is at least a{g + e). Since g < g~^, this 
quantity tends to a{g) as e \ 0. This proves for g E [p, g'^). The case 
g G {g~,p) is completely similar. Finally, when g > g^ (the case g < g~ is 
similar) we have I{g) = 00 and then Rn{g) = for large n, proving ([6]) in 
this case. □ 

Proof of Proposition \2.1\ The properties of I are clear from the definition. 

Fix ?7 G fio- In view of (jlj) and (ITOl) . the Laplace transforms of i'n{-) = 
P{{l/n)Hn = ■) have logarithmic asymptotics: 

lim - In Ep{exp{(3HniS)}) = <^{f3) - ln(2rf) 
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for all real j3. From the Gartner-Ellis theorem (Theorem 2.3.6 in [5]), the full 
statement (i) in Proposition 12.11 follows, and we obtain for open G C [0, 1] 
that 

liminf-lnz/„(G') > - inf{/(^); ^ G G H ^} , (12) 

where 

£: = {f? G [0, 1] : 3/3 Vr ^ ^, - ^\q) > (3r - ^*{r)] 

is the set of exposed points of if* from ([5]). Its complement is the set of all 
points Q such that ip* is linear in a neighborhood of q. We will improve (fT2!) 
into (ii) of Proposition 12.11 with a subadditivity argument. We start by 
showing that is differentiable at with that v^'(O) = p. Indeed, using 
Jensen inequality twice, we have 

-hiEp\e^^^"--'"P^\ < E-lnEp[e^(^"-"*')l < - InEE^pfe^^^"-"^)! 
n n n 

Computing the extreme terms and taking the limit n — > oo for the middle 
one, we get 

Q<^{(3)-(3p<\{(3)-l3p, 

which shows that v^'(O) = p since A'(0) = p. This implies that p E £ and 
that £^ is a neighborhood of p. Let g G {q~ , g'^) fl G be a non-exposed point 
of if*. For definiteness, we assume g > p. Let 

gi = sn-p{g' E £; g' < g} , g2 = mf{g' E £; g' > g} . 

Recall that if is strictly convex by Theorem 11.51 - that we will prove below 
independently. This implies that the function cannot have a linear piece 
that goes up to g^, cf. Figure [H Then, p < gi < g < g2 < g'^ , and gi, g2 & £■ 
Let 7 G (0, 1) such that g = 7f?i + (l— 7)f?2- Since the interval {gi, ^2) consists 
of non-exposed points, we have I{g) = 7/(f?i) + (1 — 7)/(f)2)- Since G is open 
and contains g, we can find e > and k,i eW such that 

, , , , ku + iv 

\u - gi\ < e,\v - ^2 < £ =^ —, ^ e G 

k + i 

with G"^ the set of r G G at distance at least e from the outside of G. The 
key fact is 
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> ^ Card G Vn{k+e) ■ ^"^^^^ ^ i^i -^,Qi+ Snk = 



X Card 1 5 G Vn(k+e) ■ Snk = x, 

Hn(k+i){S) — Hnk{S) 



ni 

> Card G Vn{k+e) ■ ^""^^^ ^ (^'i " ^'i + ^^)} 
X min Card <^ S G Vn{k+E) '■ Snk = x, 

Hn{k+t){S) — Hnk{S) 



G (^2 -e,Q2 + 



nl 



G {qi -e, 02 + e) 



Ce.Td[seVnk■■^^^eig^-s,g, + e)] 

H^:''^\s) 



X min Card |gGP„,: / ^ € (g2 - g, g2 + g)| 

||z|ji<?ifc I ni J 

with H^^'^\S) = Ylt=i vit + St + x) the Hamiltonian in the time-space 
shifted environment. Similarly, we denote by v^^i'^^ measure v^ni'^^^'^ ~ 
P{H^'''^^ G ■). The above display implies that 

n^oo n[k + Ij 

k 1 
> — — liminf — lnz/„fc((^i - e, ^1 + e)) 

k + I n^co nk 

+ — ^liminf^ min In z/^"'"'''^ ((^2 - ^'2 + 

k + i n^oo ni \\x\\i<nk ^ ' 

It is straightforward to check that 

lim inf In i/„(fc+£) (G^) < hm inf - In z/„(G) , 

n^oo n(A; + t) n^oo n 

and it is not difficult to see that 

liminf^ min \r^v"^!''''\{Q2- e.Q^^ e))> -l^Q-^) , P-a.s. (13) 

We postpone the proof of (flSll for the moment. Hence, the key inequality 
implies 

1 A; £ 

liminf-lni/„(G) > --r—l{Qx\e)-——I{Q2 + e), 
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limmf-lni/„(G) > - [7/(^?i) + (1 - 7)/(^?2)] = -/(^?), 

n^oo n 

letting e \ and k/{k+t) ^7. This yields statement (ii) in Proposition l2.1[ 
Now, let us prove ([T^ . By a standard concentration inequality (e.g., 
Theorem 4.2 in [6]), we have 

P(|lnZ„-ElnZ„| >m) <2exp| . 



Therefore we have, P-a.s. as n ^ 00, 



max 

\\x\\\<ni<n 



-lnZf'^')(/3)-^(/5) 



n 

with ziT''^^ the partition function associated to H^'^\ Since ^2 is an exposed 
point for f ll3l) follows from the Gartner-Ellis theorem. □ 

Let us comment on the above proof. We could improve (1121) into the full 
lower bound (ii) in Proposition 12. II with a subadditivity argument, implying 
convexity of the rate function. If we knew that {g^ , g'^) C £ - or, equiv- 
alently, that ip is differentiable -, we could directly conclude without this 
extra argument. We tried to prove it, but we could not. We state it as a 
conjecture: 

Conjecture 2.2 The function ip is everywhere differentiable. 

2.2 Proof of Theorem [TTS] 

1. By Jensen inequality, 

ElnZ„ < nX{p) . 

Then, f{l3) < A(/5), which implies 'p*{g) > ^*{g) from the definition of 
Legendre transform. The inequality now follows from a < —ip*. 

2. Set (PniP) = n~^ElnZ„(/5). From Theorem 1.1 in jYj we have 

for all /3 > 0. Hence, for g > p, the reciprocal functions are such that 



/5e 
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Since (A')'^ = (A*)' and {if'J-^ = «)', we have 

(A*)'(f.) < (vl'iQ) 

for all ^ > p where ip* is differentiable. Since a = —if* for ^ this 

proves the first half of the desired statement. The other half is similar. 

3. From Theorem 1.1 in |7j it is known that the set 

>V(p) = {/3 e R : ^{(5) = \m 

is an interval containing 0. Let P E W(p), and g = A'(/3) = A'(/5). From 
Theorem 2.3 (a) in 0] it is known that j3 G VV(p) implies ^p*{q) < 0. 
Then, the supremum defining A*(^) is achieved at (3, which implies the 
first equality in 

-y{Q) = -wq- xm = -ipq - m] = -v*{q) = «(f?) , 

where the second equality holds for P G W(p), the third one because 
of (p'{P) = A'(/3) = g, and the last one because (p*{g) < 0. 

Let now (3 ^ W{p), and g = X'{(3). Then, 

-X*ig) = -IPg - A(/3)] >-[Pg- > -if*{g) > aig) . 

Observe that A' is a diffeomorphism from M to (0, 1). From this we can 
identify the set V(p) defined by ([9]), 

V(p) = {A'(/3);/3g >V(p)}, (14) 

which is an interval containing p. 

4. When d = 1, it is known that W(p) = {0}, see Theorem 1.1 in [?]. 
Hence, V{p) reduces to {p}. 

5. When d > 3, from celebrated results of Imbrie and Spencer [13] . 
Bolthausen p], it is known that VV(p) contains a neighborhood of 0. In 
view of f[T^ . V(p) is in its turn a neighborhood of p. 

6. This is a consequence of [5, example 2.1.1], which shows for instance 
that, if p > TCd, then >V(p) D M^. Indeed, in view of f|T^ . this implies 
that V(p) contains [p, 1), and a is still equal to —A* at g = 1 hj upper 
semi-continuity of both functions. The case of p < 1 — tt^^ is similar. 



16 



7. This is a consequence of [51 example 2.2.1], which shows for instance 
that, if p < (l/2(i), then W(p) is bounded from above. The other case 
is similar. 

□ 

2.3 Strict convexity of the free energy 

The aim of this section is to prove Theorem II. 5[ We start with a variance 
estimate analogous to that for Gibbs field in [10]. 

Lemma 2.3 For any compact set K G there exists a positive constant 
C = Ck such that 

E{\n Zn)" {(3) >Cn , (3 e K 

Proof: The polymer measure at inverse temperature /3 with environment r) 
is the random probability measure /i„ = /i^ on the path space defined by 

/i„({5}) = exp{[3H4S)} , S = {Si, . . . , G P„ . (15) 

For simphcity we write /in(5'i, . . . , Sn) for /i„({(S'i, . . . , The polymer 

measure is Markovian (but time-inhomogeneous), and 

(lnZ„)"(/?) = Var,,„(iJ„). 

Let Si be the t-coordinate mapping on Vn given by Si(S') = St, and regard 
it as a random variable. 
Define 

I{x,y) = {z eZ'^ : \\x - z\\i = \\z -y\\i = l} , x,yeZ'^ (16) 

the set of lattice points which are next to both x and y. The set T{x, y) 
is empty except if y can be reached in two steps by the simple random 
walk from x; in this case its cardinality is equal to 2d, 2 or 1 according to 
y = X, \\y — x\\oo = 1 or \\y — x\\oo = 2. The Markov property implies 
that, under jU2n, ^1,^3, ■ ■ ■ ,^2n-i are independent conditionally on T,'^ : = 
(S2, S4, . . . , T,2n), and the law of T,2t~i given S"^ only depends on T,2t-2, ^2t, 
and has support X(J22t-2, ^2t)- 
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From the variance decomposition under conditioning, we have 



Var^,„(/72n) = E^,y^v^,SH2n I S^) + Var^,„(E^,Jif2n I S^]) 

n 

= E^,„Var^,„(5^r/(2t- l,S2t-i) | S^) 
t=i 

n 

= ^/..„Var^,„(r/(2t - 1, E^t-i) | S^) 



where £'[ • | S^], Var( ■ | S^) denote conditional expectation and conditional 
variance. To obtain the last equality we used the conditional independence. 
Define the event 

M{7],t,y,z) = {Card{r/(t,a;);a;e J(y,2;)} = 2} . 

The reason for introducing M(?7, t, y, z) is that on this event, a path 5* condi- 
tioned on St-i = y, St+i = z, has the option to pick up a rjit, St) value that 
can be either or 1, bringing therefore some amount of randomness. This 
event plays a key role here, as well as in the proof of Lemma [373] below. Note 
for further purpose that 

p(M(r^, t, y, Z))=l- (pCardlfe,.) ^ (1 _ p)CardX(,,.)) ^. _ ^) _ ^^7) 

The key observation is, for aA\ t < n and (3 & K, 

Var^,„(r/(2t-l,S2t-i) | J:') > Cl{M{r],2t - l,^2t-2,^2t)} , (18) 

where the constant C depends only on K and the dimension d. Indeed, on 
the event M{ri,2t — 1, 821-2, S2t), the variable ri{2t — l,T,2t~i) brings some 
fluctuation under the conditional law: it takes values and 1 with probability 
uniformly bounded away from provided [3 remains in the compact. Hence, 

n 

EVar^,„(ff2n) > CEY,f^2n[M{r],2t-l,^2t-2,^2t)] 



t=i 

n 



CE^Y^Y^ /U2n(S2t_2 = a:,E2t = y)l{M(r7,2t-l,x,i/)} 
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For 1 < i < n, let /in be the polymer measure in the environment 'fi{t,x) — 
'r){t, x) if t ^ i, 'rj{i, x) — for all x. Obviously, 

C-j2^\S)<US)<C+j2';:\S), SeVn, 

with positive finite C~,C~^ not depending on 71,77,/? ^ Then, with C — 

CC-, 

EVar^,„(//2n) 

n 

> C'eJ2 E &~'\^2t-2 = x,E2t = y)l{M{rj,2t - l,x,y)} 

n 

= C'eJ2 E i^?n~'\^2t-2^^,^2t^y)nM{r),2t-l,x,y)) 

n 

> 2C'p{l-p)Ej2 E £''\^2t-2^x,J:,t^y)l{\\x-y\\^<l} 

2n 

> C'p{l-p)Ej2 E ~£''\^t-2=x,J:t = y)H\\x-y\\oo<l} , 

t=2 x,y£Zd- 

since we can repeat the same procedure, but conditioning on the path at odd 
times. Finally, with ASt := St — St-i and C" = C'C^p{l — p), we have for 
alipeK,e> 0, 

2n 

EVar^,„(ff2n) > C"E^^,„Eli^^*7^^^*-i> 

t=2 

> nC"exEfX2n{An,e) , (19) 

where 

2n 
t=2 

It is easy to see that the complement 

2n 

^n,s = { E ^i^S^ = ^S^-^} > - ^)} 

t=2 
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of this set has cardinahty smaller than exp{nS{e)}, with S^e) \ as e \ 0. 
We bound 

¥{max{H2n{Sy, S G A^,J > 2nQ) < e"'^^") x Prob(i3(2n,p) > 2n^) , 

with B{2n,p) a binomial random variable. It follows that there exists some 
g{e) with g{e) \ p as e \ such that the left-hand side is less than 
exjp{—n5{eY^'^}. For all r] such that max{H2n{S)] S E ^ne} — 2?T.f?(£), we 
have the estimate 

fi2n{Ale) < exp{2n[Pg{e)-v{P) + 6{e) + o{l)]} 
< exp{2n[ip*{g{e)) + 6{e) + o{l)]} 

with o(l) ^ as n — > oo. But, as e \ 0, 

ip*igie)) + 5ie)^cp*ip) = -\ni2d)<0. 

By continuity we can choose e > such that ip*{g{e)) +6{e) < (—1/2) \n{2d), 
and E/i2n(An,e) — i> 1 as n — s> oo. Finally, from (IT^ we obtain the desired 
result for even n. The same computations apply to ^2n+ii yielding a similar 
bound. This concludes the proof of Lemma 12.31 □ 

Proof of Theorem \1.5l It follows from Lemma [2.31 that, for (3,j3' G K, 

with if'^ the right- derivative, and a similar statement for (3' < (3. Indeed, this 
inequality holds for (l/n)ElnZ„ instead of ip, and we can pass to the limit 
n ^ oc. This yields the strict convexity of (p. By a classical property of 
Legendre duality, it implies the differentiability of □ 



3 Sharp asymptotics 



Assume d > 3. Let Uq be the open set in the complex plane given by 
f/o = {/5 £ C : I Im/5| < 7r}. Then, Uq is a neighborhood of the real axis, and 
A(/3) = logE[exp{/3r7(t, a;)}] is an analytic function on Uq. Define, for n > 
and [3 G f/o. 



Wn{P) = Ep \ exp ((3 Vit, St) - nX{P) 



(20) 



t=i 
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Then, for all P G Uq, the sequence {Wn{P),n > 0) is a (^„)„-martingale with 
complex values, where Qn = (^{vitt x);t < n,x e Z^}. At the same time, for 
each n and rj, Wn{f3) is an analytic function of /5 G Uq. 
Define the real subset 



Ui = {f3eR : A(2/3) -2A(/3) < -InvTd} 



(21) 



which is an open interval (/^^ , Pi) containing (— oo < Pi < < Pi < +oo). 
The following is established in [7J: 

Proposition 3.1 Define U2 to he the connected component of the set 

|/3g?7o : A(2Re/3) -2ReA(/3) < -luTTdj 

which contains the origin. Then, U2 is a complex neighborhood of Ui. Fur- 
thermore, there exists an event Qi with P(fii) = 1 such that, 

WniP) WooiP) as n ^ 00, for all r] eQi,P , 

where the convergence is locally uniform. In particular, the limit Woo{P) is 
holomorphic in U2, and all derivatives of Wn converge locally uniformly to 
the corresponding ones of Woo. Finally, Woo{P) > for all P G Ui, F-a.s. 

For the sake of completeness we repeat the proof here. 



Proof of Proposition [X71 Since (e^) = and E[/] = E[/], we have X{P) 
\{P), and 



E 



E 



Ep[exp{pHn{S) - nX{P)}]Ep[exp{pH4S) - nX{P)}] 



= Ep02 

= EplB2 

< Ep02 

< 00 



E 



exp{pHn{S) + pHn{S) - 2nRe\{p)} 



exp I [A(2 Re/3) - 2 Re A(/3)] ^ l{St = St} 

t=i 

expUx{2ReP) -2ReX{P)]J2HSt = St} 



t=i 



(22) 
(23) 



up G U2- Indeed, the random variable Ylt^i ^{St = St} (which is the number 
of meetings between two independent (i-dimensional simple random walks) 
is geometrically distributed with parameter vr^. 



21 



For any real /? G U2, the positive martingale Wn{P) is bounded in L^, 
hence it converges almost surely and in £^-norm to a non-negative limit 
WooiP)- Moreover, the event {Woo{P) = 0} is a tail event, so it has proba- 
bility or 1. Since KWoo{P) = 1, we have necessarily Woq{P) > 0, P-a.s. 

We need a stronger convergence result. Fix a point (3 E U2 and a radius 
r > such that the closed disk D{(3,r) C t/2. Choosing R > r such that 
D{i3, R) C U2, we obtain by Cauchy's integral formula for all (3' G -D(/5, r). 



2«vr Jai5(/3 R) 2; - P' Jo 



hence 



(/3 + i?e2^™) - p' 



X^:=snp{\W,,{P')\-p' eD{P,r)}<R [ l^r.(/3 + i?e -)| ^^^ 
Letting C = {R/{R — t))^, we obtain by the Schwarz inequality 

(E[X„])2 < CE[[ |lK(/3 + i?e2'™)|2rfM] 

< Csnp{E[\Wn{P")\'];n>l,P"eDiP,R)} 

< 00 

in view of (!23l) . Notice now that X„, a supremum of positive submartingales, 
is itself a positive submartingale. Since supE[X„] < 00, X„ converges P-a.s. 
to a finite limit X^q. Finally, 

snp{\Wn{P% P' G D{p,r),n> 1} < 00 P-a.s., 

and Wn is uniformly bounded on compact subsets of U2 on a set of envi- 
ronments of full probability. On this set, {Wn,n > 0) is a normal sequence 
[2T] which has a unique limit on the real axis: since U2 is connected, the full 
sequence converges to some limit Woo, which is holomorphic on U2, and, as 
mentioned above, positive on the real axis. □ 

We do not know that Woo{P) 7^ for general P E U2, only for [5 E Ui. 
Therefore, for all 77 G fii, we fix another complex neighborhood f/3 of f/i, 
included in U2 and depending on t], such that Woo and Wn (for n large) 
belongs to C \ M_ . Recall that 

Zn{l3) = Wn{l3)e^v{n\m (24) 
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by definition. 

It is sometimes convenient to consider, for real /?, the /3-tilted law 

VnAk) = Znif^y'e^'Qnik) , G {0, 1, . . . , n} , 

which is a probability measure on the integers 0,1, ... ,n. Its mean is equal 
to (d/dp) \nZn{P), and its variance is 

D„,^ = ^lnZ„(/?). (25) 

These quantities depend also on rj, and -Dn,/3 > as soon as the Bernoulli 
configuration {r]{t,x),t < n, ||x||i < n, ||a;||i = n mod 2) is not identically 
or 1 on each "hyperplane" t = k, k = 1, . . . ,n. This happens eventually with 
probability 1, so we will not worry about degeneracy of the variance -D„,/3. 
By positivity of the variance, for all u in the range of {d/df3) lnZ„(-) there 
exists unique /3 = G M such that 

^lnZ„(/?„(n))=«. (26) 

Observe that the function /5„ is itself random. Define for /5 G M, A; G N, 

In{k) = sup{/3fc-lnZ„(/5);/3GM}-r2ln(2d) . (27) 

(We will see in the proof of Theorem 11.61 below, that /„(/c) ~ nl{k/n) with 
/ as in Proposition 12. 1[ ) For k in the range of {d/d(3) lnZ„(-), we have 

In{k) = (3n{k)k - In Zn{(3nm - nln(2d) . (28) 

Recall {13^, Pt) defined in (|2T|). 

Theorem 3.2 There exist an event with P(i72) = 1 one? a real neigh- 
borhood U4 of 0, U4 C {Pi, (3^), with the following property. Let kn G 
{0, 1, . . . ,n} he a sequence such that Pn{kn) remains in a compact subset K 
ofUi, and let Dn = Dn^p^i^k,,)- Then, for all rj G ^2, 

Qn{kn) = ^i=exp{-/„(A;„) +nln(2c/)} x (l + o(l)), 
V27rD„ 

where o(l) ^0 as n ^ 00. 
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Proof of Theorem \3.S\ Suppose that /3 is a real number. Note that the 
Fourier transform of the tilted measure is 

From the usual inversion formula for Fourier series we have 

Taking (5 = Pn{kn) and using fl28|) this becomes 

Q^(kn) = e-^"('=")+'^''^(2d) ^ _ / — ^ ^^-iknu^^ _ ^29) 

For the moment, K is any compact subset of From the Taylor 

expansion of Zn 3X (3 = (3n{kn) and (126|1 . we have 

log Z„ (^/9„(/c„) + iuj = log Zn {j3n{kn)^ + wA;„ - y^^n + Rest„ , 

where the remainder can be estimated by the Cauchy integral formula, 

|Rest„| < |u|35^=^max{| logZ„(/5')|; e D{P\5k),P" e iT} 

for all |m| < 5ki with 5/^ > equal to half of the distance from K to the 
complement of f/3. From Proposition 13.11 and the definition of f/3, the above 
maximum is less that Cxn for all n > 1, with Ck random but finite and 
independent of n. 

Moreover, in view of Proposition 13.11 and fl24ll25p . we see that 

bn = n\"{(3n{kn)) + C(/5n(A:„)) (30) 

is such that C'^n < Dn < C'l^n for some positive constants C^, C^. 

We split the integral in fl2^ according to \u\ < := (Inra/n)^/^ or not, 
and the first contribution is 



ZniPnikn) + IU 



n\<e„ Zn(^Pn{kn)^ 



du 



24 



u\<£, 



u 



exp<| - —Dn\du{l + o{l)) 



1 + 0(1)) 



27vDr 



(31) 



since SnDn'^ ^ oo by fl30|) . 

Finally, to show that the other contribution is negligible, we need the 
following fact: 

Lemma 3.3 There exist an event with P(f23) = 1, an integer random 
variable hq, a neighborhood f/5 of in M., and k > such that noijj) < 00 for 
rj E fls and 



Zn{(3 + iu) 



for Tj G ^3, [3 G f/5, u G [— vr, vr], and n > nQirf). 

With the lemma to hand, for 77, /5, u as above, we bound 



I 

J £, 



en<u<n 



Zn{Pn{kn)) 



where we have used n = 0{Dn) of fISUl) . Combined with fl^ and fl2I?l) this 
estimate yields the proof of the theorem, with fl2 = fliCifl^, and t/4 = U^nU^. 

□ 

We turn to the proof of Lemma 13.31 which states that the distribution 
i^n,/3 does not concentrate on a sublattice of Z, and is not too close from 
such a distribution. In our proof we take advantage of some (conditional) 
independance in the variables r]{t, St) under i'n,/3- This is reminiscent of a 
construction of [11] for central limit theorem and equivalence of ensembles 
for Gibbs random fields. 

Proof of Lemma \3.3[ In the notations of the proof of Lemma 12.31 



Z2n{(3) 



iuH2n 



E E 



iuH2n 
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< Eiln \E,,r, 



-^M2n -^M2n 



^iuri(2t-l,S2t-i) I 



by conditional independence of Ei, S3, . . . , T,2n-i under /i2„ given E*^. Recall 
the notation X from flTBl) and denote by 



rui = Card {x e I{S2t-2, S2t) ■ rj{2t - 1, x) = £} 



0,1, 



the number of sites which can be reached by the walk at time 2t — l and where 
?7(-) equals to and 1 respectively (mi + mo < 2(i). Then, for mo, mi > 1, 



E 



^iur]{2t+l,S2t+i) I 



mie 



l3+iu 



+ mo 



mie^ + mo 
< exp{— Cu^} , 



\u\ < n , 

where the constant C is uniform for (3 E K, and 1 < mo, mi < 2d. We obtain 



771 JuH2n 
^tJ.2n^ 



n 

^2t-2, ^2t)} \ 

t=l 



So far our arguments do not require (3 to be small. From this point, we will 
use a perturbation argument. Since iJ,2n is equal to P for /3 = 0, we study the 
term on the right-hand side for the simple random walk measure P instead of 
the polymer measure H2n, and estimate the error from this change of measure. 
This procedure is rather weak, we believe that the result of the lemma holds 
for a much larger range of (3, but we we do not know how to control the term 
in the right-hand side in a different way. 

For e > we split the last expectation according to the sum being larger 
or smaller than ne, 



E^,2^e 



iuH2n 



< e-^^"'+/i2n(5^1{M(77,2t-l,E2t_2,E20} <^ 

i=l 

n 

< e-^-' + e2"^p( ^ l{M{r], 2t-l, ^2t-2, ^2t)} < ne') (32) 

t=i 

by the obvious inequalities < H2n < 2n. For 7 G (0, 1], note that 
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Eexp[-^l{M{7],2t-l,S2t-2,S2t)}} = e-^q{S2t-2 - S,t) 

+ [1 - qiS2t-2 - S2t)] 
with q defined in f|T7|) . Then, there exists some Ci > such that 

sup [e~'^q{x) + [1 - qix)]) < exp{-Ci7} > 7 ^ (0, !]• 

x:P{T,2=x)>0, 

||x||oo<l 

Hence, 

n 

EEp exp I - 7 ^ l{M(r/, 2t - 1, E2t-2, ^2t)}} 
t=i 

n 

= Epexpj -Ci7 J]l{||S2i_2-S2i||oo < 1}} 

t=l 

= (Epexp|-Ci7l{||S2||oo<l}}^ 



2d 

with C2 > 0. Now, we choose e = C2/2, 7=1, and we get 

n 

EP(^l{M(r7,2t-l,S2t-2,S2t)} <n£) 
t=i 

n 

< e^-^'EEp exp I - 7 J] l{M(r/, 2t - 1, S2t-2, ^2*)}} 

t=i 

and then 

n 

P(P(^ l{M(r7, 2t- 1, S2t_2, ^2t)} < ne) > e""^^/^) < e""^^/^ . 
By Borel-Cantelh lemma, the set ^3 of all environments such that 

n 

P(^Y1 HM{v, 2t-l, S2t-2, ^2t)} < ne^ < e-""^^'^ eventually. 



t=\ 
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is of full measure. We define hq as the first integer (if exists) from which the 
previous bound is fulfilled, and = (— C2/4, C2/4). From (1521) we easily 
check that Lemma [3.31 holds true with k = min(C£, C2/2). □ 

Proof of Theorem \l.b\ The theorem is a corollary of Theorem 13. 2[ where VL2 
and f/3 are introduced. In particular we know that a = —rj* in t/3. Note 
that l3{g) is the maximizer in the definition of A*(^) as a Legendre transform. 
Since kn/n g, we have that /5„(fc„) — > l3{g). By fl30l) . Dn ~ nX"{P{g)), and 
by Legendre duality, 

(A*)' o A' = Id , 

and so X"{P{g)) = l/{X*)"{g). The only quantity left to be studied is In{kn)- 
Combining fl27t [24|) and performing the change of variable (3 = (3{kn/n) + v, 
we have 

In{kn) = sup{/?fc„-nA(/?)-lnl^„(/?);/5GM} 

= sup +v)kn -n\{l3{kn/n) + v) 

-hiWn{(3{knln)+v);v G ] 

\{(3{kjn)) - \{f3{kjn) + v) + \'{f3{kjn))v 

+ n\*{kn/n) 



sup \ n 

-\TiWn{(3{knln)+v);v G 
n\*{kJn)-\nWn{(3{kJn)) 



+ sup |n \{l3{kn/n)) - X{(3{kn/n) + v) + X'{(3{kn/n))v 

-\nW^{P{kJn) +v) + \nWr,{(3{kJn)y,v e 
nX*ikJn)-\nWnil3ikJn))+oil) 



= nX*{kn/n)-\nWnWQ))+o{l) 
by strict convexity of A and the fact that | ln[l^„(/5 + v)/W„ 



< \v\. 



□ 
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